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We deduce and discuss the implications of self-similarity for the stability in terms of robustness 
to failure of multiplexes, depending on interlayer degree correlations. First, we define self-similarity 
of multiplexes and we illustrate the concept in practice using the configuration model ensemble. 
Circumscribing robustness to survival of the mutually percolated state, we find a new explanation 
based on self-similarity both for the observed fragility of interconnected systems of networks and 
for their robustness to failure when interlayer degree correlations are present. Extending the self¬ 
similarity arguments, we show that interlayer degree correlations can change completely the stability 
properties of self-similar multiplexes, so that they can even recover a zero percolation threshold and a 
continuous transition in the thermodynamic limit, qualitatively exhibiting thus the ordinary stability 
attributes of noninteracting networks. We confirm these results with numerical simulations. 


I. INTRODUCTION 

Self-similarity is defined in a wide sense as the prop¬ 
erty of some systems to be, either exactly or statistically, 
similar to a part of themselves. This property is found 
in certain geometric objects that are intrinsically embed¬ 
ded in metric spaces, so that distance in the metric space 
gives a natural standard of measurement to uncover simi¬ 
lar patterns at different observation scales [T] . In complex 
networks, the definition of self-similarity is not obvious 
since many networks are not explicitly embedded in any 
physical geometry and the only available metric is the 
one induced by the collection of shortest path lengths 
between nodes. This metric has, in fact, been used to 
measure the fractal and self-similar properties of com¬ 
plex networks HIS]. However, the small-world property 
typically found in real complex networks strongly limits 
the range of scales where such properties can be observed. 

In the absence of a natural geometry, the main problem 
in the definition of self-similarity stems from the fact that 
there is, a priori, no way to decide what is the “part” of 
the system that should be compared to (and look alike) 
the “whole”. In this sense, self-similarity is not an in¬ 
trinsic property of the system but it is directly related to 
the specific procedure to identify the appropriate subsys¬ 
tem. In previous work on single networks, self-similarity 
was properly defined on the basis of a nested hierarchy 
of subgraphs and proved for general classes of models. 
These include random scale-free models with and with¬ 
out underlying metric spaces and models of growing net¬ 
works g g. Interestingly, metric network models are 
able to provide a plausible explanation for key topologi¬ 
cal properties observed in real networks ISHS], including 
scale-free degree distributions, high levels of clustering, 
the small world property, and self-similarity. 

Self-similarity has important implications in the global 
structure of networks and, in particular, in their vulner¬ 
ability to failures of their constituents. For instance, 
self-similarity alone —independently of the divergence 
of the second moment of the degree distribution— ex¬ 


plains the absence of a percolation threshold in random 
scale-free networks, with a proof that avoids the usual lo¬ 
cally tree-like and other limiting assumptions [5]. More¬ 
over, the same proof applies to ensembles of graphs with 
highly non-trivial topologies as long as they belong to 
the same self-similarity class. In g, the absence of a 
percolation threshold was also proved and numerically 
confirmed in ensembles of random networks embedded in 
metric spaces with strong clustering and in ensembles of 
growing networks with bounded topological fluctuations. 

In this work, we extend the concept of self-similarity 
to multiplexes —defined as networks of nodes intercon¬ 
nected with different classes of links, each class named 
a layer. Out of the many different self-similar ensembles 
in single networks, we chose for simplicity the configu¬ 
ration model and generalize it to multiplexes in order 
to state explicitly the definition and significance of self¬ 
similarity in such structures. In particular, we study the 
implications of self-similarity for the stability and robust¬ 
ness to failure of multiplexes with and without interlayer 
degree correlations. Circumscribing robustness to sur¬ 
vival of the mutually percolated state [9l [10], we find 
a new explanation based on self-similarity both for the 
observed fragility of uncorrelated scale-free systems of 
networks mm and for their robustness to failure when 
correlations are present [niiini. We find that interlayer 
degree correlations can change completely the stability 
of self-similar scale-free multiplexes, which can recover a 
zero percolation threshold and a continuous transition in 
the thermodynamic limit qualitatively exhibiting so the 
ordinary stability properties of single scale-free networks. 

The paper is organized as follows. In section jnj we 
review the definition of self-similarity in single-layered 
networks and extend it to multiplexes. In section m 
we discuss the self-similarity properties of the canoni¬ 
cal configuration model generalized to multiplexes, both 
with and without interlayer degree correlations. In sec¬ 
tion |IV[ we use this model to deduce and discuss the 
implications of self-similarity on mutual percolation and 
check our predictions against numerical simulations. Fi¬ 
nally, we conclude in section [V| 
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II. SELF-SIMILAR ENSEMBLES 

In the next section, we first review our findings on this 
topic in the case of single networks and, then, extend 
them to the case of multiplexes. 

A. One-layered self-similar ensembles 

Let G{{oi}) be an ensemble of sparse graphs in the 
thermodynamic limit, where {a} is the set of model pa¬ 
rameters. For example, in the case of the Erdos-Renyi 
model [TillTS] the set {a} is just the average degree (fc). 
Consider a transformation rule T that for each graph 
G e G{{oi}) selects one of G’s subgraphs. Denote the 
ensemble of these subgraphs by Gt{{c(}). The ensemble 
G{{a}) is called self-similar with respect to T if the trans¬ 
formed ensemble is the same as the original one except 
for some transformation of the model parameters, that 
is, 

GT{{a)) =G{{aT}), ( 1 ) 

where {ax} are the ensemble parameters after the filter¬ 
ing process. This definition does not assume anything 
about the transformation rule T and, in fact, the same 
ensemble can be self-similar under different rules. As 
a simple example consider the Erdos-Renyi model with 
A ^ 1 nodes and connection probability among pairs of 
nodes p = (k)/N. Now consider the transformation rule 
that selects Nx nodes uniformly at random out of the 
original N nodes, along with their connections. It is easy 
to see that such subgraph belongs to the Erdos-Renyi 
ensemble but with an average degree 

{k)T = ^{k). ( 2 ) 

Note that the average degree of subgraphs generated with 
this procedure is smaller than the average degree of the 
original network. 

In the ensembles studied in HIS] — including the stan¬ 
dard configuration model with scale-free degree distribu¬ 
tions and zero clustering, scale-free networks with finite 
clustering and metric structure, and non-equilibrium net¬ 
works, like generic growing network models—, the only 
model parameter that changes after the transformation 
is the average degree of the subgraph, {k)x- Typically, 
this average is a monotonic function of the ratio between 
the size of the original network N and the size of the 
subgraph Nx, that is. 

In this case, the sign of its derivative determines the class 
of self-similarity of the model and, in turn, the structural 
properties of the entire network. For instance, when f{x) 
is a monotonic increasing function, any graph of the en¬ 
semble contains subgraphs with an arbitrary large aver¬ 
age degree within the subgraph. This is the case of the 



FIG. 1. Illustration of a self-similar ensemble of graphs em¬ 
bedded into a metric space, a circle of radius i? ~ A, under a 
transformation that removes nodes with degrees below a cer¬ 
tain threshold [4]. In this visualization, each node is given a 
radial coordinate inversely proportional to its degree so that 
we obtain the desired subgraph by removing all nodes outside 
the blue dashed circle. 

configuration model with an scale-free degree distribu¬ 
tion with exponent 2 < 7 < 3 and, remarkably, of many 
real-world networked systems [4j. This simple property, 
together with the fact that these subgraphs belong to 
the same ensemble, imply a zero percolation threshold in 
the thermodynamic limit [5], even if 7 >> 3. Remark¬ 
ably, this is a consequence of self-similarity alone and 
not of the divergence of the second moment of the de¬ 
gree distribution. The proof in represents a powerful 
alternative to typical techniques applied to the study of 
percolation in complex networks, since it avoids the usual 
locally tree-like and other limiting assumptions. 


B. Self-similar multiplexes 

Formally, self-similarity of random multiplexes can be 
defined as for single networks. As in Eq. Q, let M{{a}) 
be a multiplex ensemble of sparse graphs in the thermo¬ 
dynamic limit, where {a} is the set of model parameters, 
now including sets of model parameters for each layer. 
Consider a transformation rule T that for each multi¬ 
plex M S M{{a}) selects one of M’s subgraphs. This 
transformation rule selects nodes in the multiplex accord¬ 
ing to specific conditions imposed on each layer. Denote 
the ensemble of subgraphs by Mx{{ct}). The ensemble 
M.{{a}) is called self-similar with respect to the trans¬ 
formation rule T if the transformed multiplex ensemble 
is the same as the original one except for some transfor¬ 
mation of the model parameters, that is, 

Mx{{a}) = M{{ax})- (4) 

To get insights on the nature and consequences of self¬ 
similarity in multiplexes, hereafter we focus on the 
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soft version of the configuration model, the simplest 
self-similar ensemble with a non-trivial degree distribu¬ 
tion [^. Nevertheless, the generalization to other ensem¬ 
bles is straightforward. 


III. THE SOFT CONFIGURATION MODEL 


The configuration model is defined as the maximally 
random ensemble of graphs with a given degree sequence, 
that is, a predefined degree assigned to each single node 
of the network [TMT5] . The soft configuration model 
(SCM) is very similar to the original one except that, in 
this case, nodes are given their expected degrees and not 
their actual degrees [niin]. This makes the model more 
appropriate to deal with structural topological correla¬ 
tions that are unavoidable when the degree distribution 
is broadly distributed [SUES]. 

In the particular case of scale-free networks, graphs 
are generated by assigning to each of the N nodes a 
hidden variable k drawn from a power-law probability 
density p{k) = (7 — 1)kq~^, k > kq. Nodes with 
expected degrees k and k' are then connected with prob¬ 
ability r{K, k') = r(/iKK'), where function r(x) < 1 is an 
arbitrary function with r(0) = 0 and r'(0) 7 ^ 0. Constant 
/i fixes the average degree {k) through the relation 


W W (5) 

Nr'{0){K.)'^ Nr'{0)KQ yy—1/ 


With this choice, it is easy to see that the average degree 
of a node with hidden variable k is proportional to k , so 
that the degree distribution scales as well as a power law 
with exponent 7 |21j . When function r(x) is chosen to 
be 


r(x) 


1 

1 -I- 1/a;’ 


( 6 ) 


the model produces maximally random graphs with a 
given expected degree sequence [23U25] . Random graphs 
with arbitrary structural correlations can be generated as 
well by choosing the appropriate connection probability 
r(x) |21] . Hereafter, we use the maximally random en¬ 
semble with connection probability given in Eq. ([^ . This 
particular ensemble has, in the thermodynamic limit, 
only two free parameters, the exponent of the degree dis¬ 
tribution 7 and the average degree (k). Notice that kq is 
a dummy parameter that can be absorbed in the defini¬ 
tion of the hidden variable k so that it can be set to unity 
at any moment. However, it is useful to keep it during 
the transformation rule that we apply below. Unlike the 
regular configuration model (where the actual degrees are 
fixed a priori) nodes in the canonical configuration model 
can end up having zero degree and, therefore, the average 
degree (k) can take any positive value, even below 1 . 

As already discussed, ensemble self-similarity is always 
tied to a particular prescription to extract subgraphs out 
of a given graph. In the case of ensembles of scale-free 


networks, the natural transformation rule selects sub¬ 
graphs by removing all nodes with degrees lower than 
a given threshold value. In the case of the SCM, the 
transformation rule T removes nodes with hidden vari¬ 
able K below an arbitrary threshold kt > ^o- In US, 
we proved that the ensemble of subgraphs so obtained 
is the same as the original one but with a transformed 
average degree 


{k)T 




(7) 


This simple result provides important insights on how 
hubs are organized within the network. We first notice 
that by varying continuously the threshold kt, we obtain 
a nested sequence of subgraphs. When 7 > 3, (fc)T is a 
monotonic decreasing function of kt- This implies that 
subgraphs made of high degree nodes are very sparsely 
connected among them. Thus, even if the original graph 
is globally connected, connectivity between two hubs is 
always mediated by chains of low degree nodes. When 
7 < 3, (fc)T is a monotonic increasing function of kt- 
In turn, this implies that, in the thermodynamic limit, 
any graph always contains subgraphs made of hubs with 
arbitrary high connectivity, even if the average degree of 
the original graph (fc) is arbitrarily small. This implies 
that such graphs always have a giant connected compo¬ 
nent and, so, the original network has a zero percolation 
threshold [5]. 


A. Generalization of the soft configuration model 
for multiplexes: Self-similarity properties 

In this paper, we restrict our analysis to self-similar 
multiplexes with two layers. Generalizations to more 
than two layers or other ensembles is again straightfor¬ 
ward. In the two-layered SCM, each node is characterized 
by two hidden variables, Ka and Kb, distributed according 
to 


P{Ka,Hb) 


1 ^ 

- p 

KaO^bO 


f Kg Kb 
\ KaO ’ KbO 


( 8 ) 


with Ka > Kao, Kb > Kbo, and p{x,y)dxdy = 1. 

In this way, (Ka) and (Kb) are proportional to parame¬ 
ters Kao and Kbo so that they can be set to unity at any 
moment. In each layer, pairs of nodes connect with con¬ 
nection probabilities ra{p,aKaKa) and rb{pbKbK'j^), where 
parameters pa and pb read 

_ (^a) 1 _ i^b) 

" Nr'MiKa)^ “ Nr',{0){Kb)^- 

Notice that the only relation between the two layers 
comes from the joint distribution p{Ka,Kb), which may 
encode interlayer degree-correlations. 

As for the transformation rule T, analogously to the 
case of single networks, given a multiplex generated from 
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this ensemble, we remove nodes in the multiplex such 
that their hidden variables Kq and Kb in each layer are 
below certain threshold values KaT and Kbx- Next, we 
analyze under which conditions the multiplex SCM is 
self-similar. 


1. Self-similar scale-free multiplexes with uncorrelated 
interlayer degrees 


When Ka and Kb are uncorrelated variables, the joint 
degree distribution corresponds to the factorization of 
the degree distributions of each layer, so that self-similar 
ensembles of subgraphs can only be achieved if the one- 
layer degree distributions are scale-free, that is, 

P(a;,2/) = (7a - l)(7b - (10) 


Thus p{x, y) is the factorization of two homogeneous 
functions of degrees — 7 a and — 76 , which gives a bi- 
dimensional homogeneous function of degree —a = 
— {"fa + 7b)- After the transformation, the remaining 
nodes in the subgraph are distributed according to the 
same scale-free distributions once we replace Kao —>■ KaT 
and Kbo —> KbT- The number of nodes that remain in the 
subgraph is 


Nt 



KbO 

KbT 


76-1 


N. 


( 11 ) 


The transformation does not change neither the hidden 
variables of filtered nodes nor their connection probabil¬ 
ity, which implies that parameters pa and Pb remain in¬ 
variant in the subgraph. Therefore, by combining Eq. ([^ 
and Eq. (11), we conclude that the transformed ensemble 


is self-similar with re-scaled average degrees 


and 


{ka)T = 


{kb)T = 


' KbO \ 1 

( KaT 

KbT/ 

^KaO 

KaO V“~^ 1 

( KbT 

KaT J 

V KbO , 


3-7a 


(ka) ( 12 ) 


3-76 


(kb). (13) 


Notice that in multiplexes with uncorrelated degrees the 
two thresholds, KaT and KbT, are completely indepen¬ 
dent. 


2. Self-similar scale-free multiplexes with correlated degrees 

In multiplexes with correlated degrees, self-similarity 
is achieved when the joint distribution p{x, y) is a bi- 
dimensional homogeneous function of degree a, that is, 

p{ax,ay) = a~°‘p{x,y) Vo. (14) 

When the degrees in each layer are correlated, this con¬ 
dition enforces a relation between the two thresholds. 


i. e. KaxlKao = KbT/Kbo, which are not independent 
anymore [26) . Using the homogeneity property Eq. (141, 


it is easy to check that the number of nodes within a 
subgraph with Ka > KaT and simultaneously Kb > KbT = 
KbOKaT/KaO IS 


Nt 


- [ Y 

KKaT/ 


N. 


(15) 


As in the case of uncorrelated multiplexes, the transfor¬ 
mation does not change neither the hidden variables of 
filtered nodes nor their connection probability, which im¬ 
plies that parameters pa and pb remain invariant in the 
subgraph. Then, by combining Eq. ([^ and Eq. (151 we 
conclude that the ensemble is self-similar with re-scaled 
average degrees in each layer 


(fco)r = 


KaT 

^aO 


-| 4—a 


(ka) and {kb)T = 


KaT 


“I 4—O' 


KaO 


{kb). 

(16) 


IV. STABILITY OF MUTUALLY PERCOLATED 
STATES IN SELF-SIMILAR SCALE-FREE 
MULTIPLEXES 

As mentioned in the introduction, the stability prop¬ 
erties of systems of networks can be radically different as 
compared to single networks depending on the patterns 
of connectivity between layers [siiniiii]. We shall show 
that self-similarity can explain several of the previous re¬ 
sults on the robustness of systems of networks and can 
predict new behaviors in a large class of self-similar mul¬ 
tiplexes. Notice that the results presented here are qual¬ 
itatively valid in multiplex ensembles beyond the SCM if 
those present similar self-similarity properties. 

We study stability in terms of the robustness of the per¬ 
colated state. In multiplexes, the percolated state can be 
defined according to different criteria. Here, we assume 
that nodes in each layer mutually depend on nodes in 
other layers and that only the nodes that belong to the 
giant mutually connected component remain functional. 
The giant mutually connected component of a multiplex 
network (MCC) is defined as the largest set of nodes that 
are mutually connected by at least one path in each layer 
traversing nodes in the MCC lidn]- 

For single networks, perturbations in the form of a ran¬ 
dom failure of a fraction of 1—p nodes produce typically a 
critical phase transition for a specific value Pc, so that be¬ 
low Pa the network is fragmented into small components. 
In multiplexes with a MCC, perturbations can propagate 
back and forth between the layers so that even small ini¬ 
tial failures can produce avalanches of damage leading to 
a discontinuous collapse of the MCC [5]. Site percolation 
on random multiplexes has shown indeed a discontinu¬ 
ous hybrid transition at some finite value of the num¬ 
ber of nodes removed, where the size of the MCC drops 
abruptly to zero, like in a first order transition, while the 
critical behavior is only observed above the transition. 
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like in a second-order one PE]. So, perturbations are 
amplified by the interaction between the layers and sys¬ 
tems of networks are said to be more fragile as compared 
to single networks. The presence of interlayer degree cor¬ 
relations can however revert the situation m- Interde¬ 
pendent networks with mutually dependent nodes hav¬ 
ing identical degrees are statistically more robust than 
random coupled networks with the same degree distri¬ 
bution. Besides, when 7 < 3, they disintegrate via a 
second-order phase transition —in the same way as non¬ 
interacting networks— and are thus very resilient against 
random failures m- More structured systems of corre¬ 
lated interconnected networks or with overlaps have been 
proved to be robust to failure as well [131 ES] . 

Next, we assess the resilience of MCCs to random fail¬ 
ures in scale-free multiplexes on the basis of their self¬ 
similarity properties and check numerically our predic¬ 
tions. Before that, we note that the average degree {k) 
in the SCM ensemble defined in Sec. |TTT] is equivalent to 
the site percolation probability p and it can then be used 
in robustness studies as the control parameter. Indeed, 
when a random fraction of 1 — p nodes is removed from 
a given graph of the ensemble, the hidden variables ks 
of the remaining nodes are distributed as in the original 
graph and the connection probability among them re¬ 
mains unchanged. However, the number of nodes in the 
subgraph is pN. Since p remains unchanged, Eq. ([^ 
implies that this ensemble is self-similar under a ran¬ 
dom removal of nodes with a modified average degree 
{k)T = p{k). This means that, in the thermodynamic 
limit, removing a random fraction of nodes 1 — p of a 
network with average degree (k) is equivalent to gener¬ 
ating a graph of the same ensemble but with an average 
degree p{k). Because of this equivalence, hereafter we use 
{k) as the control parameter of the percolation properties 
of the ensemble. 



A. Fragility of uncorrelated scale-free multiplexes 
explained by self-similarity 


Single scale-free self-similar networks in the thermody¬ 
namic limit with 7 < 3 always contain subgraphs made of 
hubs with arbitrary high connectivity, even if the average 
degree of the original graph (fc) is arbitrarily small, which 
implies that such graphs always have a giant connected 
component and, so, a zero percolation threshold |S] . This 
makes such structures robust to random failures. In the 
case of uncorrelated multiplexes, the question is whether 
it is still possible to find a continuous set of nested sub¬ 
graphs such that the average degrees within the sub¬ 
graphs increase in both layers simultaneously. In that 
case the multiplex would be robust to random failures, 
being able to maintain a MCC despite perturbations. 

To have a nested ensemble of subgraphs, HbT must 
be either constant or a monotonic increasing function 
of KaT (or vice versa). Let = giKar) be such func¬ 
tion. Then, the condition for Eqs. (12) and (13) to be 


FIG. 2. Sketch of a self-similar multiplex with two layers, A 
and B, with uncorrelated degrees across layers. In panel a, 
At is A’s subgraph made of nodes with high degrees, such 
that the giant component of layer A is almost contained in 
At. The set of nodes in At induces a subgraph in layer B, 
Bt. In panel b, because of the self-similarity properties of 
the ensemble, Bt is similar to B but, due to the impossibil¬ 
ity to satisfy simultaneously the inequalities in Eq. (171, it 
has a smaller average degree than B. Thus, its giant compo¬ 
nent —which is the potential candidate to be the MCC of the 
multiplex— is also reduced. In panel c, this process can be 
iterated at infinitum and can lead to the fragmentation of the 
multiplex. 


simultaneously monotonic increasing functions of KaT is 


7a - 1 ^ KaTg'jKaT) ^ 3 - 7 a 
3 - 7b KaOgiKar) 7b “ 1 ' 


(17) 


However, these inequalities can only hold if the lower 
bound is smaller than the upper bound, which is equiva- 
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lent to the inequality a = 7a + 7b <4. This is clearly not 
possible in scale-free sparse graphs with 7 a and 76 in the 
range (2,3), implying that, while it is possible to have 
a sequence of subgraphs with increasing average degree 
in one of the layers (if one of the inequalities is satis¬ 
fied), the same sequence of subgraphs has necessarily a 
decreasing average degree in the other layer. 

This result explains the fragility of scale-free systems 
of networks first reported in [ 3 . In single scale-free 
networks, global connectivity is mainly provided by the 
interconnection of high degree nodes, which gives the 
main explanation for their robustness. In uncorrelated 
scale-free multiplexes, the situation is different. Our self¬ 
similarity argument starts by selecting a subgraph of high 
degree nodes in layer A and so an almost fully connected 
subgraph that contains the majority of nodes of the giant 
component of layer A (see panel a in Fig.[^. However, as 
our previous result shows, the average degree in layer B 
of the subgraph induced by the subgraph in A is smaller 
than in B and, thus, its giant component in B —which 
is the candidate set to contain the MCC of the mutually 
percolated multiplex— is also reduced. We could now 
select a subgraph of the subgraph in layer B such that 
its average degree is high enough to contain its layer B 
giant component. However, the average degree of the 
induced sub-subgraph in layer A will decrease below its 
original value, and so its giant component. This pro¬ 
cess can be iterated at infinitum and, at each iteration, 
the size of the potential subgraph to contain a MCC is 
reduced. We thus conclude that the MCC cannot be sus¬ 
tained by high degree nodes alone and must rely on the 
connectivity of low degree nodes. This makes scale-free 
multiplexes always more fragile than more homogeneous 
networks with the same average degree. 


B. Robustness of correlated scale-free multiplexes 
explained by self-similarity 


The picture changes completely when the degrees in 
each layer are positively correlated. In the case of sparse 
scale-free self-similar multiplexes with uncorrelated de¬ 
grees in the two layers, a = ja + lb > 4 so that the 
conditions for a stable MCC are not fulfilled. However, 
when Ka and Kb are positively correlated, it is possible to 
find ensembles with 3 < a < 4. As an example, consider 
the joint distribution 




7(7-1)2^-^ 

{x + yp+^ 


(18) 


Its marginal distribution is p{x) = (7 — 1)2'*'“^(1 -|- 
x)~'^ [29]. From here, the conditional average is {x\y) = 
{y + l)/{l — 1 ), so that the correlation between x and y 
increases when 7 —>■ 2. The joint distribution Eq. (18) 
is a homogeneous function with a = 7 -I- 1. Therefore, 
according to Eqs. ( |I^ , when 7 < 3 the ensemble has 
self-similar subgraphs with increasing average degree in 
both layers simultaneously. This readily implies that the 



FIG. 3. Average degrees (left column) and the size of the 
largest connected components (right column) as a function of 
filtering parameters KaT and Kbx- Panels a and b show re¬ 
sults for a multiplex network with uncorrelated degrees where 
KaT = KbT- Panels c and d show results for a multiplex net¬ 
work with uncorrelated degrees where k^t is fixed to the mini¬ 
mum value Kbo . Panels e and f show the results for a multiplex 
network with correlated degrees where KaT = KbT- Solid l ines 
correspond to the analytical results given by Eqs. (12 I, (131 
and (I 6 I. In all cases, the multiplex network is composed of 
two layers with A = 5 x 10® nodes, 7 = 2.8, (kmin) = 2, and 
we evaluated the absolute size of the largest connected com¬ 
ponents SaT and SbT in individual layers A and B, the size of 
the MCC St, and the size of the network Nt after applying 
the corresponding transformation. 


ensemble always possesses a MCC so that its percolation 
threshold is zero in the thermodynamic limit. Besides, 
the “transition” is continuous, in the sense that the rela¬ 
tive size of the MCC approaches zero monotonously when 
p —?► 0. This generalizes the result found in [T^ for net¬ 
works with identical degrees in both layers and makes an 
important step forward as it quantifies the precise level 
of correlations (and so the value of a) that is needed to 
go from a hybrid discontinuous transition to a continuous 
one. 
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C. Numerical simulations 


To check numerically the predicted stability proper¬ 
ties of self-similar scale-free multiplexes, we generated 
two-layered multiplexes using the canonical configura¬ 
tion model. In all cases, N = 5 x 10® and {kmin) = 2. 
For uncorrelated scale-free multiplexes, we used the joint 
probability distribution Eq. (10), while we implemented 
correlations according to Eq. (18). Finally, to compute 
mutually connected components, we implemented an ef¬ 
ficient algorithm based on Ref. EO], which keeps track of 
all the MCCs, not only the giant, present in a multiplex. 
The algorithm represents each layer of the multiplex by 
the dynamic connectivity structure defined in |31j . This 
structure allows for maintaining information about net¬ 
work components and their sizes, while updating a graph 
by deletion or insertion of edges. The algorithm works in 
two phases. First, we find MCCs in the initial multiplex 
and second, we calculate the size of the giant MCC for 
all values of the parameter p. 

To compute all MCCs in the initial multiplex, we iden¬ 
tify connected components for each layer separately and 
if needed, we reconnect all single components by adding 
a minimum number of ad hoc edges. Thus, after this 
step every layer is a single connected component and the 
multiplex a single MCC. Next, we sequentially delete all 
add hoc edges. Each single removal creates two separated 
components in the given layer. We then check all possible 
node pairs, where each node in the pair belongs to a dif¬ 
ferent component and remove, in all other layers, edges 
connecting them. Whenever any removed edge breaks 
a connected component into two, we have to continue 
with the removal of all edges that connect disconnected 
components in all other layers. Finally, when all ad hoc 
edges are removed, all layers consist of connected compo¬ 
nents corresponding to MCCs. In the second phase, we 
generate a random sequence defining the order of node 
removals. Removal of each node is accomplished by re¬ 
moving all its adjacent edges from all layers. Every edge 
is removed in the same way as ad hoc edges in the first 
phase of the algorithm. Similarly as in the first phase, 
after removing the node all layers consist of connected 
components corresponding to MCCs. The size of the 
largest component is outputted as the size of the largest 
MCC for the corresponding p value. 


In Fig. we show the average degrees in the sub¬ 
graphs and the size of the largest connected components 
in each layer and the MCC as a function of the filter¬ 
ing thresholds KoT and In all cases, networks are 
scale-free with 7 = 2.8. In uncorrelated multiplexes, the 
average degrees of the subgraphs cannot increase simul¬ 
taneously as the thresholds increase. This is shown in 
Fig. [3] a for KaT = i^bT and in Fig. ic for KbT = = 1 - 

As clearly seen in the figures, the only possibilities are 
that the average degrees decrease simultaneously (when 
KaT = kiit) or that the average degree of one of the lay¬ 
ers increases while the other decreases (when k^t is con¬ 
stant). This induces the fragility of the MCC which, as 



FIG. 4. Comparison between the percolation properties of 
scale-free multiplexes with uncorrelated (left panel) and cor¬ 
related (right panel) degrees. Panels a and b show the rel¬ 
ative size of the largest mutually connected component vs 
the fraction p of nodes remaining undamaged. In both cases 
A = 5 X 10® and 7 = 2.2, 2.5, and 2.8. Each curve corre¬ 
sponds to one complete random sequence of node removals. 
Panels c and d show the susceptibility x as a function of site 
occupation probability p for scale free multiplexes of 7 = 2.8 
and different sizes. The different curves x(p) are computed 
from 10'* complete random sequences of node removals. In all 
cases, multiplexes are composed of two layers and {kmin) = 2 . 


shown in Fig.j^b and d, reduces its size abruptly at some 
relatively small value of the threshold. Interlayer degree 
correlations change completely the picture. In Fig. e 
and f, we show the average degrees in the subgraphs and 
the size of the different components for KaT = ^br in a 
canonical configuration model multiplex ensemble with 
the joint degree distribution given by Eq. (18). In this 
case, it is possible to produce sequences of subgraphs with 
increasing average degrees in both layers simultaneously, 
so that the MCC becomes very robust. Finally, the inset 
in Fig. f shows the relative size of the MCC (relative 
to the remaining number of nodes after the filtering pro¬ 
cess), which approaches I for large values of the thresh¬ 
olds, indicating that, as predicted, such self-similar mul¬ 
tiplex contains a small but macroscopic subgraph that is 
completely connected in both layers simultaneously. 

To get further insights into the percolation proper¬ 
ties of self-similar multiplexes, we adopt the conventional 
percolation criterion of measuring the breakdown of the 
largest MCC. We computed the relative size of the largest 
MCC versus the fraction of nodes p remaining in the mul¬ 
tiplex for different values of the power-law exponent 7 . 
Results are shown in Fig.j^a for multiplexes with uncor¬ 
related degrees and in Fig. |^b for correlated ones. For 
all values of 7 , the transition between the mutually per¬ 
colated and the fragmented states is discontinuous in the 
uncorrelated case while it is continuous and approaching 
zero in the correlated case. This can be corroborated by 















the scaling of the susceptibility vs the system size, where 
the susceptibility x is defined as 


(S^) - (S)^ 

(S) 


(19) 


Here S is the size of the largest MCC at any value of 
p and averages are taken over a large number of com¬ 
plete random sequences of node removals. This quantity 
is able to distinguish between discontinuous, continuous, 
and hybrid phase transitions. In continuous phase tran¬ 
sitions, X shows a clear peak close to the critical point 
that diverges as the system size increases. Instead, in 
discontinuous transitions, x shows a discontinuity at the 
critical point but no dependence on the system size. In 
the case of hybrid phase transitions, x shows a diverg¬ 
ing peak approaching the critical point from one side, a 
discontinuity and then a size independent behavior on 
the other side. According to these criteria, Fig. [i] c in¬ 
dicates that the transition is hybrid in multiplexes with 
uncorrelated degrees whereas Fig. d indicates that y 
has a continuous divergence with a peak that approaches 
zero in the thermodynamic limit. This clearly corrobo¬ 
rates our theoretical prediction about a zero percolation 
threshold but with critical fluctuations when p —> 0. 


whose average degrees grow with their depth in the hier¬ 
archy [5] and without the need of usual limiting require¬ 
ments. 

In this paper, we have extended the concept to mul¬ 
tiplexes and illustrated its importance by assessing the 
stability of scale-free multiplexes in terms of their self¬ 
similarity properties. To state in a clear and explicit 
way the definition and relevance of self-similarity, we 
have focused on the SCM ensemble. However, we should 
stress that the results presented here are qualitatively 
valid in other multiplex ensembles with similar features, 
that is, with similar self-similarity properties, degree dis¬ 
tributions and interlayer degree correlations. Interest¬ 
ingly, the observed fragility of scale-free multiplexes or 
the robustness to failure of correlated systems of networks 
can be explained and predicted based only on their self¬ 
similarity characteristics. In particular, we have found 
that scale-free multiplexes can recover a zero percolation 
threshold and a continuous transition in the thermody¬ 
namic limit, and so the ordinary stability properties of 
single scale-free networks. Self-similarity can as well have 
important implications for other critical phenomena tak¬ 
ing place in multiplex structures when the critical point 
is a function of the connectivity of the system. 


V. CONCLUSIONS 

Self-similarity is a widespread property in network 
models and has also been observed in many real world 
networks |3]. Beyond the mathematical beauty of self¬ 
similarity, this property has important implications for 
the structural properties of networks. The power of the 
concept was illustrated in single-layered networks by the 
proof of a zero percolation threshold for a general class 
of self-similar networks, which only required the self¬ 
similarity property with a hierarchy of nested subgraphs 
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